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Abstract. We have investigated the build up of anisotropic flows within a parton cascade approach 
at fixed shear viscosity to entropy density Tj /s to study the generation of collective flows in ultra- 
relativistic heavy ion collisions. We present a study of the impact of a temperature dependent r\/s{T) 
on the generation of the elliptic flow at both RHIC and LHC. Finally we show that the transport 
approach, thanks to its wide validity range, is able to describe naturally the rise - fall and saturation 
of the V2{pr) observed at LHC. 
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INTRODUCTION 

The RHIC program at BNL has shown that the azimuthal asymmetry in momentum 
space, namely the elliptic flow V2 , is the largest ever seen in HIC suggesting that an 
almost perfect fluid with a very small shear viscosity to entropy density ratio, rj/s, 
has been created [1]. The first measurement at LHC in Pb+Pb at 2.767W [2] shows 
that the integrated elliptic flow as a function of collision energy increase of about 30% 
compared to the flow measured at RHIC energy of 200 GeV, while the V2(pr) measured 
at LHC comparted to that of RHIC does not change indicating an increase in the average 
transverse momentum. It remains to be understood if this means an equal rj/s of the 
formed plasma or it is the result of different initial conditions and possible larger non- 
equilibrium effects. 

The most common approach to study viscous correction is viscous hydrodynamics 
at second order in gradient expansion according to the Israel-Stewart theory [3, 4, 5]. 
This approach has been implemented to simulate the RHIC collisions providing an 
upper bound for rj/s < 0.4. Such an approach, apart from the limitation to 2+1D 
simulations, has the more fundamental problem of a limited range of validity in r\/s 
and in the transverse momentum pj. In these proceedings we discuss results within 
the relativistic transport approach that has the advantage to be a 3+ ID approach not 
based on a gradient expansion in viscosity that is valid also for large rj/s and for 
out of equilibrium momentum distribution allowing a reliable description also of the 
intermediate pj range. In this pj region viscous hydrodynamics breaks its validity 
because the relative deviation of the equilibrium distribution function 8f/f eq increases 
with pj becoming large already at pr > 3T ~ IGeV . In the following we will show 
the results obtained with a parton cascade approach where the EoS is fixed to the 
one of a free massless gass e — 3P = and the mean free path A is finite. A more 



quantitative comparison with the experimental data would require the inclusion of mean 
field dynamics associate to an equation of state P(e) according to 1QCD results [6]. 
A first step in this direction has been discussed in [7, 8] while the implementation of 
a quasi-particle model [9] in the contest of transport theory is in progress along lines 
similar to [10]. 



THE PARTON CASCADE AT FIXED r\ /s 

Our approach is a 3 + 1 dimensional Montecarlo cascade [11] for on-shell partons 
based on the stochastic interpretation of the transition rate discussed in Ref. [12]. In 
kinetic theory under ultra-relativistic conditions the shear viscosity can be expressed 
as r\ = (4/15)p < p > X with p the parton density, X = [pc tr ]~ l the mean free path 
and < p > the average momentum. Therefore considering that the entropy density for a 
massless gas is s = p(4 — fi/T), /i being the fugacity, we get: 

where o tr is the transport cross section. In our approach we solve the relativistic Boltz- 
mann equation with the constraint that r\/s is fixed during the dynamics of the collisions 
in a way similar to [13] but with an exact local implementation a more detailed discus- 
sion of the method is in [11]. In fact fixing r\ / s we can evaluate locally in space and 
time the strength of the cross section o tr (p,T) needed to have r\/s at the wanted value 
by mean of the following formula: 

4 <p> 1 
tr 15 p (4-^/7)7] /s K) 

This approach is equivalent to have a total cross section of the form Ojot = 
K(p,T)o p qcd > GpQCD where K takes into account the non perturbative effects 
responsible for that value of viscosity. Note that this approach have been shown to 
recover the viscous hydrodynamics evolution of the bulk system [4, 11], but implicitly 
assume that also high pj particles collide with largely nonperturbative cross section. 
We show here that both at RHIC and LHC there are signatures of the disappearence of 
the large non perturbative physics with increasing pj. 



Effect of temperature dependent ri/s(T) 

In our calculation the initial condition are longitudinal boost invariant with the initial 
parton density dN/dr\(b = 0) = 1250 at RHIC and dN/dr\(b = 0) = 2250 at LHC. 
The partons are initially distributed in coordinate space according to the Glauber model 
while in the momentum space at RHIC (LHC) the partons with pj < po = 2GeV 
(Pt < Po = 4GeV) are distributed according to a thermalized spectrum with a maximum 
temperature in the center of the fireball of 2Tq (3.57c), while for p T > po we take the 
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FIGURE 1. Different temperature dependent parametrizations for r]/s. The orange area take into 
account the quasi-particle model predictions for 77 /s [9]. 



spectrum of non-quenched minijets according to standard NLO-pQCD calculations. We 
also start our simulation at the time to = 0.6/m/c at RHIC and to = 0.3/m/c at LHC. 

In order to study the effect of the kinetic freezeout on the generation of the elliptic 
flow we have performed two calculations one with a constant Anr\ /s = 1 during all 
the evolution of the system (red dashed line of Fig.l) the other (shown by black solid 
line in Fig.l) with Anrj/s = 1 in the QGP phase and an increasing rj/s in the cross over 
region towards the estimated value for hadronic matter Aitr\ /s = 8 [ 14] . Such an increase 
allows for a smooth realistic realization of the kinetic freeze-out. In Fig. 2 it is shown the 
elliptic flow V2(pt) at mid rapidity for 20% — 30% centrality for both RHIC Au+Au at 
= 200GeV (left panel) and LHC Pb+Pb at y/s = 2.767W (rhigh panel). As we can 
see at RHIC energies, left panel of Fig. 2, the V2 is sensitive to the hadronic phase and 
the effect of the freeze out is to reduce the V2 of about of 25%, from red dashed line to 
black solid line in left panel of Fig. 2. For the pj range shown we get a good agreement 
with the experimental data for a minimal viscosity r}/s ~ 1/ (An) once the f.o. condition 
is included. At LHC energies, right panel of Fig. 2, the scenario is different, we have that 
the V2 is less sensitive to the increase of 77/5 at low temperature in the hadronic phase. 
The effect of large r\/s in the hadronic phase is to reduce the V2 by less than 5% in the 
low pj region, from red dashed line to the black solid line in right panel of Fig. 2. This 
different behaviour of V2 between RHIC and LHC energies can be explained looking at 
the life time of the fireball. In fact at RHIC energies the life time of the fireball is smaller 
than that at LHC energies, 5/m/c at RHIC against the about 10 fm/c at LHC. Therefore 
at RHIC the elliptic flow has not enough time to fully develop in the QGP phase. While 
at LHC we have that the V2 can develop almost completely because the fireball spend 
more time in the QGP phase. 

Due to this large life time of the fireball at LHC and the larger initial temperature is 
interesting to study the effect of a temperature dependence in f] / s. In the QGP phase 
rj/s is expected to have a minimum of rj/s ~ (An)~ l close to Tq as suggested by 1QCD 
calculation [15]. While at high temperature quasi-particle models seems to suggest 
a temperature dependence of the form rj/s ~ T a with a ~ 1 — 1.5 [9]. To analyze 
these possible scenarios for r\ /s in the QGP phase we have considered two different 
situation one with a linear dependence Anrj/s = T/Tq = (e/eo) 1//4 (blue line) and the 
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FIGURE 2. Differential elliptic flow V2(pr) at mid rapidity for 20% — 30% collision centrality. On the 
left panel, the orange band indicate RHIC results measured by STAR and the orange points on the right 
panel are the LHC results measured by the ALICE collaboration, data taken by [2]. The red dashed line is 
the calculation with 4^tj /s = 1 during all the evolution of the fireball and without the freeze out condition, 
while the black blue and green lines are calculations with the inclusion of the kinetic freeze out and with 
Akt] /s= 1, 4-717] /s °< T and 4-7tr] /s °< T 2 in the QGP phase respectively. 

other one with a quadratic dependence 4^7] /s = (T /To) 2 = (e/£o) 1//2 (green line) where 
£0=1 -7GeV I fm is the energy density at the beginning of the cross over regions where 
the rj/s has its minimum, see Fig.l. At RHIC energies the V2 is essentially not sensitive 
to the dependence of r\ js on temperature in the QGP phase, see the blu and green lines 
in the left panel of Fig. 2. However the effect on average is to decrease the value of V2 but 
at low px < 1.5 GeV the vi{pj) appears to be insensitive to 7]/s(T) while a quite mild 
dependence appears at higher pj where however the transport approach tends always to 
overpredicted the elliptic flow observed experimentally. At LHC energies the build-up 
of V2 is more affected by the r}/s in the QGP phase and on average it is reduced of about 
a 20%. In any case still a strong temperature dependence in r]/s has a small effect on 
the generation of V2 we found that with a constant or at most linearly dependent on T 
rj/s the transport approach can describe the data at both RHIC and LHC at least up to 
pj ~ 2GeV. However the transport approach should keeps its validity also at higher px, 
but as previously said, the agreement with data seem to weaken at pj > 2GeV both at 
RHIC and LHC. We discuss the underlying reasons in the next section. 



Impact of high p T partons on v 2 

In our approach we have that Ojot = Ko p qcd therefore we have large cross section 
independently of the pj of the colliding particles. But we know that particles with 
high energies should collide with the pQCD cross section. In order to take into account 
the proper scattering cross section for hard collisions we extend our previous approach 
allowing for a K factor that depends on the invariant energy of the collision K{s) which 
gives the connection between the non pertubative interacting bulk and the asymptotic 
pQCD limit. We choose this function in such a way that at high energies K(s) — >■ 1 and 
we get the correct assumed pQCD limit. For the function K(s) we choose an exponential 





FIGURE 3. Left: v 2 (pr) at mid rapidity for 20% - 30% collision centrality at RHIC. The dashed lines 
are the calculations with K = const and for 4nr] /s = 1 and Anr\/s °< T with f.o. respectively for black 
and blue curves while the solid lines are the same but with K(s/A 2 ). Right: V2{pr) at mid rapidity and for 
20% — 30% collision centrality at LHC with the same legend, data taken from [2]. 



form K(s/A 2 ) = 1 + ye~ s l^ , where A is a scale parameter that fix the energy scale 
at which the pQCD behaviour begins to be reached. While y plays the same role of 
K in the old calculations and it is determined again in order to keep fixed locally the 
rj/s. Therefore we can repeat the same procedure as described in the previous section 
but now with <3j ot = K(s/A 2 )o p qcd- Due to its physical meaning we of course expect 
A to be greater than 2GeV, in particular we have performed different calculation for 
different value of A and we have obtained that for A > AGeV the elliptic flow becomes 
less sensitive to the value of the parameter A. Specifically in our calculation we have 
considered the value A = AGeV . As we can see at RHIC energies, left panel of Fig. 3, we 
have that K(s/A 2 ) does not affect at all the V2(pt) f° r Pt < 2GeV, in other words high 
pr parton at RHIC energies does not affect the generation of the V2 of the bulk. Instead 
we have a reduction of the V2 for pj > 3GeV and with the inclusion of K(s/ A 2 ) the V2 
becomes a decreasing function of pj for pj > 3GeV in perfect agreement with what is 
observed experimentally (orange band). In Fig. 3 (right) we compare in a large range 
the V2(pr) at LHC energy with (solid) and without (dashed) the inclusion of an energy 
dependent K factor and for two T dependence of the f] / s. We notice that the two sets of 
experimental data refer to different method of V2 measurements, namely V2 [2] (circle) 
and V2[4] (square) and our theorethical results should be compared to V2[4] because 
event-by-event fluctuations are not considered. As we can see at LHC energies the V2 
is sensitive to K(s/A 2 ) already at pr ~ l.SGeV quite lower than the RHIC case, in other 
words the many high pj partons that we have at LHC energies affect the generation of 
the V2 of the bulk. Similar results we have when we include a 7]/s(T) in the QGP phase. 

At low pr the raise of the V2 is an effect of a strong interacting fluid with a very small 
viscosity. In this regime we have that particles with low pj interact non perturbatively 
with large cross sections and therefore we get a description in agreement with hydrody- 
namics. With the increase of the pr of the partons the pQCD limit begins to be important 
and for pj > 3 — AGeV the elliptic flows starts to be a decreasing function of pj. The dis- 
appearance of the non perturbative effect significantly affects the V2(pr) making faster 



and stronger (~ 20 — 25%) the fall in the elliptic flow in the range 3GeV < pj < SGeV . 
Finally for pj > SGeV in our calculation seems to appear the saturation of the simi- 
larly to the experimental data and typical of a path-lengh mechanism as in jet quenching 
models [16]. An analysis with better statistics is required. In this range of pj the only 
effect is that given by the pQCD limit. 

CONCLUSION 

We have investigated within a transport approach at fixed r\ /s the effect of a temperature 
dependent r\/s at RHIC and LHC energies. At RHIC we have seen that the elliptic 
flow is more sensitive to the kinetic freeze out (hadronic phase) and still of 25% of 
V2 depends on it. At LHC we get an opposite effect, nearly all the V2 comes from the 
QGP phase and the r\/s of the hadronic phase is irrilevant. We get for both at RHIC 
and LHC a good agreement with the data when the ratio r\/s 1 /(4k) or a linear T 
dependence is considered, in general we observe not a large sensitivity of V2 to the T 
dependence in rj /s. Furthermore we have seen that at LHC the large ammount of particle 
with p T > 4GeV/c interacting nearly perturbatively cannot be neglected. The interplay 
between perturbative and non-perturbative behaviour seems to have an important effect 
on the generation of V2 at intermediate pj and it could explain the rapid raise and fall of 
vi{pt) in < pr < SGeV/c shown in the experiments. 
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